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Abstract 

We present a comparison of the strong couplings of light {u, d, and s), 
c, and b quarks determined from multijet rates in flavor-tagged samples of 
hadronic decays recorded with the SLC Large Detector at the SLAC Linear 
Collider. Flavor separation on the basis of lifetime and decay multiplicity 
differences among hadrons containing light, c, and b quarks was made using 
the SLD precision tracking system. We flnd: af'^/af = 0.987 ±0.027(stat) ± 
0.022(syst) ± 0.022(theory), a^Jaf = 1.012 ± 0.104 ± 0.102 ± 0.096, and 
al/af = 1.026 ± 0.041 ± 0.041 ± 0.030. 
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A fundamental assumption of the theory of strong interactions, Quantum 
Chromodynamics (QCD), is that the strong couphng is independent of quark flavor. 
This can be tested by measuring the strong coupling in events of the type e+e~ (l(i{g) 
for specific quark flavors q. Although an absolute determination of as for each quark flavor 
would have large theoretical uncertainties it is possible to test the flavor-independence 
of QCD precisely by measuring ratios of couplings in which most experimental errors and 
theoretical uncertainties are expected to cancel. Since it has recently been suggested that 
a flavor-dependent anomalous quark chromomagnetic moment could modify the probability 
for the radiation of gluons, comparison of the strong coupling for different quark flavors may 
also provide information on physics beyond the Standard Model. 

Comparisons of for 6 or c quarks with for all flavors made at PETRA [Q were limited 
in precision to ±0.41 (c) and ±0.57 (6) due to small data samples and limited heavy quark 
tagging capability. LEP measurements of a^/a"'^*'^ have reached precisions between ±0.06 
and ±0.02 under the assumption that is independent of flavor for all the non-6 quarks. 
The OPAL Collaboration has measured a{ / af^ for all flve flavors / with no assumption on 
the relative value of for different flavors |^ to precisions of ±0.026 for h and ±0.09 to 
±0.20 for the other flavors. The kinematic signatures used to tag c and light quarks suffer 
from low efficiency and strong biases, due to preferential tagging of events without hard 
gluon radiation. 

The SLC Large Detector (SLD) at the SLAC Linear Collider (SLC) is an ideal 
environment in which to test the flavor independence of strong interactions, tracking 
capability of the Central Drift Chamber (CDC) and the precision CCD Vertex Detector 
(VXD) P], combined with the stable, micron-sized beam interaction point (IP), allows 
us to select bb{g) and — > qiqi{g) {qi = u, d, s) events using their quark decay 

lifetime signatures with high efficiency and purity and low bias. Here we present the flrst 
precise measurements of a^/af^, a^/af^, and a^'^'^/af^ using this technique, and making no 
assumptions about the relative values of a^, and a^''^. 
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This analysis is based on the 1.8 pb~^ of e^e~ annihilation data collected during the 
1993 run of the SLD at the SLC at a mean center-of-mass energy of ^/s = 91.26 GeV. The 
trigger and selection criteria for hadronic decays are described in Ref . |1|] . The efficiency 
for selecting a well-contained QQig) event was estimated to be above 96% independent 

of quark flavor. The estimated 0.10 ± 0.05% background was dominated by Z^ t^t~ 
events. This analysis used charged tracks measured in the CDC and in the VXD |[l|. 

We used normalized impact parameters d/cTd as the basis for quark flavor tags, where d 
is the signed distance of closest approach of a charged track to the IP in the {x-y) plane 
transverse to the beam axis, and 0"^ is the error on d. A resolution on d of 10.8 /im has been 
measured using fJ'^f^' decays, and the spatial resolution on the average transverse IP 

position has been measured to be 7 /im 0. The distributions of d and d/aa are modeled 
well by the SLD simulation 0]. Tracks used for event flavor tagging were required to have: 
at least one VXD hit; at least 40 CDC hits, with the first hit at a radius less than 39 cm; 
a combined CDC+VXD fit quality y/2x^ — V2^d.o.f. — 1 < 8.0; momentum greater than 
0.5 GeV/c; < 250 ^m; and to miss the IP by less than 0.3 cm in the x-y plane and by less 
than 1.5 cm in z. Tracks from candidate and A decays and 7-conversions were removed. 

Figure 1 shows the distribution of nsig, the number of tagging tracks per event with 
d/cTd > 3. The data are well described by a Monte Carlo simulation of hadronic Z^ decays 
||T0| with parameter values tuned ||lT| to hadronic e~^e~ annihilation data, combined with 
a simulation of the SLD. For the simulation, the contributions of events of different quark 
flavors are shown separately. The leftmost bin contains predominantly events containing 
primary u, d, or s quarks, while the rightmost bins contain a pure sample of events containing 
primary b quarks. The event sample was divided accordingly into three parts: those events 
with risig = were deflned to be the ucis-tagged sample; those with 1 < risig < 3 were the 
c-tagged sample; and those with rigig > 4 were the 6-tagged sample. The efficiencies e for 
selecting events (after cuts) of type i {i = uds, c, b) with tag i, and the fractions 11 of events 
of type i in the i-tagged sample, were calculated from the Monte Carlo simulation to be: 
{e,U)uds = (77%, 86%); {e, U), = (59%, 38%); and {e, n)^ = (46%, 94%). 
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Figure 1. The measured distribution of the number of tracks per event with 
d/(Td > 3.0 (points). The histograms show the flavor composition estimated from 
a Monte Carlo simulation (see text). 

Jets were then reconstructed using iterative clustering algorithms. We used the 'E', 'EO', 
'P', and 'PO' variations of the Jade algorithm, as well as the 'Durham' ('D') and 'Geneva' 



('G') algorithms |T^. We divided events into two categories: those containing: (1) two jets, 
and (2) three or more jets. The fraction of the event sample in category 2 was defined as 
the 3-jet rate R^. This quantity is infrared- and colinear-safe and has been calculated to 
0{af) in perturbative QCD |[T^JT3[] . For each algorithm, the jet resolution parameter Uc 
was chosen to be as small as possible subject to the requirement that 0{a1) QCD provides 
a good description of -R3 measured in our global sample of all flavors p|,p^. This choice 
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maximizes while avoiding the 'Sudakov region' at low i/c where multiple gluon emission 
requires that large logarithmic terms of 1/yc be resummed in order to describe the data [|1| . 
The resulting yc values are listed in Table 1. 
Table 1. Results for R'^/Rf^, derived from Eq. 1; see text. Errors shown are statistical. 



Algorithm 


Vc 




R^/Rf 


Rl/Rf 


factor Rl/RjiactoT 


E 


0.080 


0.941±0.042 


1.212±0.173 


0.980±0.063 


0.995 


0.958 


EO 


0.050 


0.975±0.036 


1.113±0.145 


0.981±0.053 


0.994 


0.945 


P 


0.030 


1.001±0.027 


0.985±0.109 


1.007±0.041 


0.992 


0.929 


PO 


0.030 


1.014±0.026 


0.899±0.102 


1.037±0.039 


0.992 


0.929 


D 


0.015 


0.989±0.035 


1.096±0.145 


0.947±0.049 


0.991 


0.921 


G 


0.030 


1.032±0.020 


0.942±0.079 


0.952±0.030 


0.989 


0.915 



The -R3 for each of the j quark types [j = uds, c, b) was extracted from a maximum 
likelihood fit to and ng, the number of 2-jet and 3-jet events, respectively, in the z-tagged 
sample: 

^3 = E (^|^3)^^3 + ^1^3) (1 - ^^3)) fN . (1) 

Here is the total number of selected events corrected for the event selection efficiency, 
and is the Standard Model fractional hadronic width for decays to quark type j. The 
matrices £^(2-+2) ^(3~»3) efficiencies for an event of type j, with 2- or 3-jets at the 

parton level, to pass all cuts and be tagged as a 2- or 3-jet event, respectively, of type i. 
Matrices £(2_»3) and £^('^_^2) efficiencies for an event of type j, with 2- or 3-jets at the 

parton level, to pass all cuts and be tagged as a 3- or 2-jet event, respectively, of type i. 
This formalism explicitly accounts for modifications of the parton-level 3-jet rate due to 
hadronization, detector effects, and tagging bias. These matrices were calculated from the 
Monte Carlo simulation. The efficiencies for correctly tagging a 2-jet event and a 3-jet event 
differ by an average of 5.7%, 8.3%, and 30.3% for the uds, c, and b tags, respectively. 
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Equations 1 were solved using 3-jet events defined by each of tlie six algoritfims. Tlie 
ratios R\/ Rf^-, wfiere Rf^ is tlie 3-jet rate in the total event sample, are shown in Table 1. Av- 
eraged over all six algorithms the correlation coefficients from the fit are: uds-c : —0.76, uds- 
h : 0.30, c-h : —0.55. The statistical errors were calculated using the full covariance matrix. 

The 3-jet rate in heavy quark (6, c) events is expected to be reduced relative to that in 
light quark events by the diminished phase-space for gluon emission due to the quark masses. 
We evaluated the suppression factors, Ry R^ and R'^/R'^, for each jet algorithm and i/c value 
according to Ref. [|l5l, assuming b (c) quark masses of 4.75 GeV/c^ (1.50 GeV/c^). These 
factors are listed in Table 1, and were used to correct the measured 3-jet rate ratios. 

To 0{al) in perturbative QCD, Rsivc) = ^(l/c)"^ + {B{yc) +C{yc))a1, where the C(a^) 
coefficient includes a term B{yc) from 3-parton states calculated at next-to- leading order, 
and a term C{y^ from 4-parton states calculated at leading order. Hence, the ratio of the 
strong coupling of quark type j to the mean coupling in the sample of all flavors, a{/af^, 
can be determined from: 

Rfiyc) A{y,)af + [B{y,) + C{y,)] {af^ ' ^ ^ 

where A{yc), B{yc), and C{yc) for the different jet-finding algorithms were evaluated using 
Refs. ll^Jl^. Using our measured values of af^(M|) determined from jet rates [Q, we found 
that for the E, EO, P, PO, and D algorithms, the leading-order QCD calculation C{yc)(y'^ lies 
below the experimental 4-jet rate by roughly a factor of two. We increased C{yc) ad hoc for 
these algorithms, so as to describe the data. Equation 2 was solved to obtain ai/af^ for 
each jet algorithm; the results are shown in Fig. 2. The errors include contributions from the 
statistical error, as well as the experimental systematic errors and theoretical uncertainties. 

We considered systematic effects that could modify the tagging efficiencies. In each case 
the error was evaluated by varying the appropriate parameter in the Monte Carlo simulation, 
recalculating the matrices e, performing a new fit to Eq. 1, and rederiving ai/af^. Suitable 
variation about the world average value of each parameter was considered [H. The errors 
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Figure 2. Values of ai/af^ derived for each of the six jet algorithms for each of the 
quark flavors j (see text). The error bars on the averages include the statistical and 
systematic errors and the total theoretical uncertainty. 

are summarized in Table 2, where averages over the six algorithms are shown. The largest 
contributions result from limited knowledge of the heavy quark fragmentation functions and 
B decay multiplicity. The uncertainty in BR[Z^ cc) also produces large variations in 
o-l/af^ and /af^. Contributions from h hadron lifetimes, the fraction of in B meson 
decays, h baryon production rates, and the charm hadron decay multiplicity are small. The 
detector systematic error is dominated by the uncertainty in the charged track reconstruction 
efficiency. No systematic variation of the results was found when the event selection cuts, 
tag criteria, or Hc values were changed. 
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Table 2. Contributions to the systematic error on ai / af . 



Son rpp 


A ( a""^" 1 a^"""^ 


A f / a^"^ 




b physics 


0.008 


0.060 


0.033 


c physics 


0.017 


0.060 


0.011 


Detector modehng 


0.003 


0.032 


0.017 


Monte Carlo statistics 


0.011 


0.048 


0.014 


QCD uncertainty 


0.003 


0.011 


0.012 



We considered sources of uncertainty in the QCD predictions that affect the values of 
ai/af^ derived from Eq. 2. For each jet algorithm these include variation of the QCD 
renormalization scale within the range allowed by our measurements of jet rates in the global 
sample []T1| and variation of the heavy quark masses used in the phase-space correction 



factors by ±0.25 GeV/c^. In addition, the shifts in a{/af^ due to the ad hoc increase of 
the coefficient C{yc) were conservatively assigned as an uncertainty. The variation of the 
results due to uncertainties in parton production and hadronization were found to be small. 
These contributions were added in quadrature to yield the total QCD uncertainties listed in 
Table 2. 

There is significant scatter among the ai/af^ values derived from the different jet 
algorithms. In order to quote a single ai/af^ value for each flavor j, we made the conservative 
assumption that the results are completely correlated, and we calculated the unweighted 
mean values and errors over all six algorithms. We obtained 

^ = 0.987 ± 0.027 (stat) ± 0.022 (syst) ± 0.022 (theory) , 

^ = 1.012 ± 0.104 (stat) ± 0.102 (syst) ± 0.096 (theory) , 

^ = 1.026 ± 0.041 (stat) ± 0.041 (syst) ± 0.030 (theory) , (3) 
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where the theoretical uncertainty is the sum in quadrature of the QCD uncertainty from 
Table 2 and the rms of the results over the six algorithms. These averages are also shown 
in Fig. 2. The variation of results among jet algorithms, presumably due to different 
uncalculated 0{al) QCD contributions, dominates the theoretical uncertainty, is not small 
compared with experimental errors, and has not been considered in previous analyses [^Q. 

In conclusion, we have used hadron lifetime information to separate hadronic decays 
into three flavor samples with high efficiency and purity, and small bias against events 
containing hard gluon radiation. From a comparison of the rates of multijet events in these 
samples, we find that the strong coupling is independent of quark fiavor within our sensitivity. 
These are the first such results using a precision vertex detector for fiavor separation at 
the Z^. This represents the most precise test for uds events. Our findings are consistent 
with measurements performed at LEP using different flavor-tagging techniques [^,|^ . 

We thank the personnel of the SLAC accelerator department and the technical staffs of 
our collaborating institutions for their outstanding efforts. We thank C. Ng, T. Rizzo, and 
E. Maina for their helpful contributions. 
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